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Abstract 

Quasinormal modes have played a prominent role in the discussion of pertur- 
bations of black holes, and the question arises whether they are as significant as 
normal modes are for self adjoint systems, such as harmonic oscillators. They can 
be significant in two ways: Individual modes may dominate the time evolution of 
some perturbation, and a whole set of them could be used to completely describe 
this time evolution. It is known that quasinormal modes of black holes have the first 
property, but not the second. It has recently been suggested that a discontinuity in 
the underlying system would make the corresponding set of quasinormal modes com- 
plete. We therefore turn the Regge-Wheeler potential, which describes perturbations 
of Schwarzschild black holes, into a series of step potentials, hoping to obtain a set of 
quasinormal modes which shows both of the above properties. This hope proves to 
be futile, though: The resulting set of modes appears to be complete, but it does not 
contain any individual mode any more which is directly obvious in the time evolution 
of initial data. Even worse: The quasinormal frequencies obtained in this way seem 
to be extremely sensitive to very small changes in the underlying potential. The ques- 
tion arises whether - and how - it is possible to make any definite statements about 
the significance of quasinormal modes of black holes at all, and whether it could be 
possible to obtain a set of quasinormal modes with the desired properties in another 
way. 

PACS codes: 97.60.Lf, 04.30.-fx, ll.lO.Qr, 02.60.+y 
Paper: |Kr-qc/9602032 



I. Introduction 

Quasinormal modes are single frequency modes dominating the time evolution of per- 
turbations of systems which are subject to damping, either by internal dissipation or by 
radiating away energy. Due to the damping, the frequency of a quasinormal mode must 
be complex, its imaginary part being inversely proportional to the typical damping time. 
Examples are oscillations of stars damped by internal friction. In general relativity, damp- 
ing occurs even without friction, since energy may be radiated away towards infinity by 
gravitational waves. This leads to a new class of quasinormal modes of neutron stars which 
are not present in the Newtonian case [1,2]. Recently it has also been suggested [3,4] to 
describe the behavior of leaky optical cavities using a quasinormal mode formalism, rather 
than studying the surrounding infinite system that the cavity is embedded in. 
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Even linearized perturbations of black holes exhibit quasinormal modes [5], despite 
the absence of any oscillating fluid. The characteristic modes of black holes have first been 
found by numerical calculations [6] , leading to attempts to develop a formalism to describe 
characteristic oscillations without having to refer to a specific initial perturbation, based 
only on the properties of the underlying system. 

The study of linearized perturbations of Schwarzschild black holes leads to the Schro- 
dinger-like equation (e.g. [7]) 

^|;"{x) + (u^ - V{x)^ ^|;{x) = 0, (1) 

where ifj represents the wave function of the perturbation and V{x) is the so-called Regge- 
Wheeler potential 

I originates from the expansion of the perturbation in terms of spherical harmonics, and r 
is considered a function of x such that x = r + ln(r — 1). All units are chosen such that 
G = 1 and c = 1, the radial coordinate is scaled such that the horizon is at r = 1, i.e. 
2M = 1, where M is the mass of the black hole. 

Quasinormal modes are formally defined by those solutions of Eq. (1) which satisfy 
purely outgoing boundary conditions, i.e. 



iIj{x) ""^^ e+^'^^ and ^l){x) ""^^ e"*^^. (3) 

We have chosen the time dependence to be e*'^*, quasinormal frequencies of stable systems 
must therefore have a positive imaginary part. 

In the case of a Schwarzschild black hole, it is known [8] that there is an infinite number 
of discrete frequencies to which allow solutions of Eqs. (1) - (3). However, the significance of 
these frequencies and their associated modes is not immediately clear. In the case of normal 
modes, the underlying mathematical problem is self adjoint, and the normal modes form a 
complete set of solutions. Therefore, the behavior of the system can be analyzed completely 
in terms of its normal modes. There is no such theorem for quasinormal modes, since Eqs. 
(1) - (3) do not define a self adjoint problem. The completeness of a set of normal modes 
also makes it possible to define unambiguously what we mean by the excitation strength of 
a given normal mode as a result of some initial data. Again, this definition does not simply 
carry over to quasinormal modes, not even to the fundamental ones which are obvious in 
the time evolved data. 

Nevertheless, quasinormal modes are believed to play a significant role for the study of 
perturbations of black holes. This view is supported by the fact that the fundamental quasi- 
normal frequency (the one with the smallest imaginary part) is generally found to dominate 
the time evolution of some initial perturbation. Even a few of the more strongly damped 
modes, which are not immediately obvious in the time evolved data, can be extracted with 
suitable techniques [9] . 
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However, the fundamental quasinormal mode can only be an approximation of the 
actual time evolved perturbation. An interesting question is whether the set of all quasi- 
normal modes of a black hole can be used, in analogy to a set of normal modes of a self 
adjoint system, to completely describe the behavior of a perturbation of the black hole, 
i.e. the time evolution of some initial data for the perturbation equation. Quasinormal 
modes cannot be complete in the usual sense, meaning that at any given time, the solution 
can be represented, over all space, as a sum of quasinormal modes. Due to the boundary 
conditions (3), solutions corresponding to damped oscillations have to grow exponentially 
as a; — ±00. 

However, it may be possible to represent the time dependence of a solution at a given 
point in space completely as a sum over quasinormal modes. A model systems with this 
property has been studied by Price and Husain [10]. In general, however, quasinormal 
modes do not form such a complete system. This is related to Eq. (1), together with the 
boundary conditions (3), not defining a selfadjoint problem. Therefore, the question of the 
mathematical as well as physical meaning of quasinormal modes and of their usefulness in 
practical computations has to be studied for every system in particular. 

Ching et al. [11] have argued that a system which is described by a Schrodinger-like 
equation such as Eq. (1) will indeed have a complete set of quasinormal modes if the 
potential V{x) has a discontinuity either in the potential itself or in any of its derivatives. 

On the other hand, it is known that the quasinormal modes of a Schwarzschild black 
hole (or of a neutron star) cannot be complete. One reason is that the late time behavior 
of perturbations of a black hole is first dominated by damped oscillations ("quasinormal 
ringing"), but at very late times it exhibits a power-law tail [9]. Such a tail cannot be 
obtained by a sum of quasinormal modes. In addition, the set of quasinormal frequencies 
of black holes does not contain frequencies with large real parts, i.e. oscillations with very 
short periods. Therefore, any initial data corresponding to such oscillations cannot be 
adequately represented by quasinormal modes. 

The argument of Ching et al. [11] therefore implies that if the (smooth) Reggc- Wheeler 
potential is changed slightly such that a discontinuity, however small, is introduced in the 
potential or any of its derivatives, the set of quasinormal modes of the black hole will 
change from incomplete to complete. On the one hand, this would be very convenient, 
since it would allow the study of arbitrary perturbations in terms of this expanded set of 
quasinormal modes, just as the time dependent behavior of any self adjoint system can 
be studied in terms of its normal modes. On the other hand, one would not expect the 
physical response of the black hole, i.e. the time evolution of some perturbation, to be 
significantly affected by a small change in the potential. We are therefore facing a possible 
contradiction: The actual time evolution of some initial data will probably not be affected 
significantly by this small change in the potential, while the quasinormal mode spectrum, 
which is generally believed to represent crucial aspects of this time evolution, should change 
drastically. 

In an attempt to shed some light on this apparent contradiction, we will replace the 
Reggc- Wheeler potential in Eq. (1) by a step potential. This is a potential which is 
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piecewise constant, constructed in such a way that it approximates the original, continuous 
Regge- Wheeler potential. Of course, we expect any given step potential to have a quasinor- 
mal modes spectrum which is somewhat different from that of the smooth Regge- Wheeler 
potential. Therefore, we construct a whole series of such step potentials and study the 
behavior of their quasinormal frequencies as the number of steps is increased, i.e. as the 
step potentials are allowed to approximate the continuous Regge- Wheeler potential better 
and better. 

II. Procedure 

In the following, we will always use the Regge-Wheeler potential (2) with / = 2. There 
are, of course, many different ways to construct step potentials which eventually approxi- 
mate a given smooth potential. We will present the results for three such possibilities: 

(i) The difference in potential between steps is (roughly) constant. In order to achieve 
this, the length of the steps is variable and depends essentially on the derivative of the 
smooth potential. 

(ii) The length of the steps is constant, while the difference in height between the steps 
depends on the smooth potential. However, the length of the steps may be chosen 
differently to the left and to the right of the maximum of the potential if the potential 
is not symmetric around its maximum. 

(iii) As in (i), but an exponential damping is applied to the smooth potential before the 
steps are constructed. 

In all cases, the value of the step potential for a particular step is generally determined by 



^step(a^n-i <x<Xn) = / V{x)dx for n = -Nst + 1 . . . A^st, (4) 

where A^st is the number of steps to either the left or the right of the maximum of the 
potential. For method (i), the position XNst of the last jump on the right is determined by 



^step(a^Arst-i < < XNst) = 777^ = / V(x)dx, (5) 



2A^St XNst - XNst-l 



with a corresponding condition for the last step on the left. This ensures that 



/+00 r+00 
Vstepix)dx= / Vix)dx. (6) 
-00 J — 00 

This condition is not satisfied by methods (ii) and (iii). 

Figure 1 (i) - (iii) shows the resulting step potentials for different numbers of steps. 
Method (i) obviously leads to the last step on the right of the potential's maximum be- 
coming extremely long. This is due to the slow decay of the Regge-Wheeler potential for 
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a; — > oo: V{x) as x ^ oo, while V{x) ~ exp(a;) as x 



— oo, as can be verified 



from Eq. (2). The last step being so much longer compared to all other steps seems rather 
artificial. In fact, with increasing A'^st, the last steps grows longer, while the steps near the 
maximum of the potential become shorter, as one would indeed expect for a refinement of 
the step potential. 

Therefore, in method (iii) we have imposed an exponential damping on the potential 
before constructing the steps: 



Note that this leads to an exponential decrease with the same power as for x — oo. 
Of course, a^cut itself has to grow with increasing A^st if the step potential is supposed to 
approximate the smooth potential arbitrarily well for Nq^ oo. We use Xcut '■— ^^max + 
4(iVst/4)V4^ where the maximum of the Regge- Wheeler potential is at Xmax — 1.1947. 

A related problem occurs in method (ii): Here, the last step would become much higher 
than the differences in height of the preceding steps, again due to the slow decrease of the 
Regge- Wheeler potential for x ^ oo. In order to avoid this, several steps before the last 
are reduced in height such that the step height approaches linearly. Again, this artificial 
decrease sets in at larger values of x if the number of steps increases. 

The determination of the quasinormal modes of these step potentials is quite simple: 
On the nth step, the solution is given by il;{xn-i < x < Xn) = A^e*'^"^ + S„e~*'^"^, with 
kn = yJiJ^ — Vstep(a^)) and determined according to the boundary conditions. For 

a given (complex) frequency a;, we start with V'+oo = e"*'*'^ to the right of the potential 
and, moving to the left, determine and at every step by the standard procedure of 
matching i\){x) and '(/''(x) at every step. If — if^-ooix) = e^'^^^ after reaching the left 
side of the support of the potential, i.e. if B-^^^ = 0, the frequency a; is a quasinormal 
frequency. 



The quasinormal frequencies determined in this way are shown in Fig. 2. The details 
of the behavior of the step potentials' quasinormal frequencies, such as the value of the 
fundamental frequency, the imaginary parts of the 'lined up' frequencies, and even their 
spacing, depend on the method of construction. The dominating features, however, are 
independent of it. In all cases, the spacing between the frequencies becomes closer as the 
number of steps increases, but they generally occupy the same part of the complex plane. 
In particular, they remain 'lined up' more or less parallel to the real axis, rather than to 
the imaginary axis as the quasinormal frequencies of the smooth Regge- Wheeler potential. 
There is no indication that the frequencies cease to reach arbitrarily large real parts, i.e. 
arbitrarily small oscillation periods, as the number of steps becomes very large. 

As A'^st increases, there are individual frequencies 'escaping' from the line of frequencies 
towards the imaginary axis. However, even these frequencies do not seem to approach the 
quasinormal frequencies of the Regge- Wheeler potential. 
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III. Results 
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Figure 2: Quasinormal frequencies for step potentials constructed using methods (i) - (Hi). 
The quasinormal frequencies of the smooth Regge- Wheeler potential are included for com- 
parison. 
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(i) Variable step size 


(ii) Fixed step 


size 


(iii) Exponential cutoff 


smooth 


Index 


4 steps 


64 


256 


4 steps 


64 


1024 


4 steps 


64 
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1 


0.329 


0.084 
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0.756 
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0.006 
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0.292 


0.104 


0.030 


0.232 


0.259 


0.123 




4 


0.727 


0.202 


0.102 


1.101 


0.574 


0.156 


1.079 


1.111 
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0.298 


0.131 
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2.007 


1.351 






0.180 
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0.419 
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0.043 


0.468 


0.380 


0.284 




20 


2.673 
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0.427 
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2.060 


0.635 
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3.502 


2.296 






0.231 


0.013 


0.021 


0.550 


0.271 


0.051 


0.599 


0.131 


0.310 




50 


6.460 
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1.517 
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5.125 
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0.091 
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0.710 
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0.093 


0.737 


0.527 


0.377 
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2.027 
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ol 
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0.749 


1.610 
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0.951 
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0.295 
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0.385 






0.411 




0.432 


0.442 


0.957 


o4 




3.991 
0.236 


0.391 
0.457 






0.496 
0.473 






0.335 
0.471 




o5 






3.997 
0.229 

















Table 1: Quasinormal frequencies of step potentials for all three methods of construc- 
tion and various step sizes. Each entry lists the real part of the frequency as the upper left 
and the imaginary part as the lower right number. The first three frequencies of the smooth 
Regge- Wheeler potential are included for comparison. 



Table 1 lists some quasinormal frequencies for all three methods of construction and 
different numbers of steps, along with the 'outlying' frequencies and the fundamental quasi- 
normal frequencies of the smooth Regge- Wheeler potential. In the case of method (i), the 
long step at large x causes numerical instabilities. Therefore, results for A'^st = 1024 are 
not available; the frequencies for Nst = 256 are listed instead. 

In order to check whether the behavior we found could be the result of some numerical 
artifact, we have explicitly calculated the time evolution of some Cauchy data subject to 
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the time dependent version of Eq. (1) 

inserting step potentials constructed using method (i), with A^st = 1, -^St = 4, and A^st = 
64, for V{x). The initial data consisted of a bell shaped peak with compact support, 
incident onto the potential from the left, the point of observation is at a; = 120. The result 
for A'^st = 4 is shown in Fig. 3, and it manages to combine two seemingly contradictory 
requirements: It is only slightly different from the time evolution data for the continuous 
Regge- Wheeler potential (included in Fig. 3 for comparison), and it also represents the 
totally different spectrum of quasinormal frequencies of the step potential. The similarity 
is seen shortly after the initial burst has passed (Fig. 3(i)); actually, the time evolution 
here is dominated by the fundamental quasinormal frequency of the smooth potential, 
which is not present in the spectrum of the step potential. At very late times, however, the 
perturbation oscillates with the fundamental frequency of the step potential (Fig. 3(ii)), 
while the data for the smooth potential is damped more strongly and is therefore invisible 
in Fig. 3(ii). However, this difference concerns only a tiny fraction of the total energy 
radiated by the system. 





Stationary calculation 


Time evolution 




1 step 


64 steps 


1 step 


64 steps 


Index 


Re(c<;) Im(a;) 


Re(a;) Im(a;) 


Re(a;) Im(a;) 


Re(a;) Im(a;) 


1 
2 
3 
4 


0.5836 0.0165 
0.6786 0.0550 
0.8187 0.0975 
0.9859 0.1351 


0.0841 0.0064 
0.1203 0.0151 


0.5828 0.0169 
0.6792 0.0566 
0.8245 0.1022 
1.0333 0.1374 


0.0841 0.0064 
0.1203 0.0152 



Table 2: Comparison of quasinormul frequencies obtained from the stationary calcula- 
tion described in this paper with frequencies obtained by a fit of time evolution data. 



We determined the dominating frequencies of these late time oscillations by fitting 
quasinormal modes to the time dependent solution, their frequencies not being fixed, but 
parameters of the fit procedure. The fundamental frequencies obtained in this way are given 
in Tab. 2, they agree very well with the fundamental quasinormal frequencies obtained from 
the stationary calculation. Some deviation is expected, since the time dependent calculation 
works on a fixed grid and does not know about the exact location of the jumps in the step 
potential, or even about the potential being a step potential at all. 

IV. Discussion 

Introducing discontinuities in the Regge- Wheeler potential has a significant influence 
on the quasinormal mode spectrum, even if the jumps become very small. Most of the 
frequencies become 'lined up' roughly parallel to the real axis, as they do for the simple 
square barrier potential, rather than being lined up along the imaginary axis like the 
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Figure 3: Time evolution of initial data, observed at x — 120; using the smooth Regge- 
Wheeler potential and the step potential with 4 steps and construction principle (i). 

quasinormal frequencies of the smooth Regge- Wheeler potential. This is a necessary, but 
not sufficient, condition for the quasinormal modes forming a complete set for the solutions 
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of the perturbation equation. 

Time evolution computations confirm that the very late time behavior of some initial 
perturbation is indeed dominated by the fundamental quasinormal mode of the step poten- 
tial under consideration, rather than by the fundamental mode or the tail belonging to the 
smooth potential. However, due to the exponential damping, this very late time behavior 
represents only a tiny fraction of the total response of a black hole, or some other system, 
to an initial perturbation. Most of the total response occurs at earlier times, where it is 
dominated by quasinormal ringing at the fundamental frequency of the smooth potential. 

Our results suggest that it is indeed possible to obtain an expanded set of quasinormal 
modes which can completely represent the time evolution of perturbations of black holes. 
This conclusion increases the significance of quasinormal modes for the study of perturba- 
tions of black holes. On the other hand, we also lose something: None of the frequencies of 
the new modes is obvious in the largest part of the time evolution of the perturbation. It 
is generally assumed that at least the fundamental quasinormal frequencies have physical 
meaning in the sense that they will dominate the time evolution of a perturbation of the 
black hole, as it is known for the quasinormal frequencies of the continuous Regge- Wheeler 
potential. The quasinormal spectrum of a very similar step potential, on the other hand, 
might actually lead to a complete system of quasinormal modes, but there is no single mode 
or frequency which has an obvious relationship to the time dependence of the perturbation. 

We have also seen that even a very small deviation from the original potential leads 
to a completely different quasinormal mode spectrum. It appears rather farfetched to 
assign significance to quantities which are so sensitive to small details, rather than to 
the overall features of the problem. On the other hand, this sensitivity might well be 
a consequence of some specific property of the step potentials, as it does not seem to 
occur for other cases that have been studied: Ferrari and Mashhoon [12], for example, 
have replaced the Regge- Wheeler potential by a Poeschl-Teller potential, approximating 
the Regge- Wheeler potential mainly around its maximum. They achieve good agreement 
between the fundamental quasinormal frequencies of the two potentials. 

How, then, can we know if the quasinormal spectrum of any given system does indeed 
tell us something about the time evolution of a perturbation of this systems? What are the 
criteria which distinguish a "meaningful" from a "meaningless" spectrum? Are quasinormal 
frequencies always strongly sensitive to small changes in some parameters of the problem? 
And can we find another way to obtain an "expanded" set of quasinormal modes which is 
complete, but leaves the "original" modes more or less intact? 

There is an intriguing parallel between the quasinormal spectrum we have found for 
step versions of the Regge- Wheeler potential and those of neutron stars: The frequencies 
of the so-called wave modes, which are associated mainly with the metric outside the star, 
rather than with the fluid inside, show a picture rather similar to what we have found 
here for the quasinormal frequencies of the step potentials, with the frequencies of the w- 
modes corresponding to the frequencies lined up along the real axis, and the frequencies 
of the wii-modes corresponding to the 'outlying' frequencies [2]. A better understanding 
of neutron star perturbations as well as Schrodinger or wave equations with discontinuities 
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is required to decide whether this paraUel is a coincidence or not. It is possible, however, 
that the junction conditions at the surface of the neutron star provide an analogy for some 
discontinuity in a potential. 

The present work shows that the significance of the quasinormal mode spectrum of a 
black hole is not yet clearly understood. The most intriguing questions may be whether it 
is possible to "complete" the set of quasinormal modes in such a way that the fundamental 
modes still represent the major features of the time evolution, and how the excitation of a 
particular quasinormal mode can be measured, or even be given a unanimous meaning. 
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